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Abstract 

The nature of the evolution of the magnetic field, and of current sheet formation, at three- 
dimensional (3D) magnetic null points is investigated. A kinematic example is presented which 
demonstrates that for certain evolutions of a 3D null (specifically those for which the ratios of 
the null point eigenvalues are time-dependent) there is no possible choice of boundary conditions 
which renders the evolution of the field at the null ideal. Resistive MHD simulations are described 
which demonstrate that such evolutions are generic. A 3D null is subjected to boundary driving 
by shearing motions, and it is shown that a current sheet localised at the null is formed. The 
qualitative and quantitative properties of the current sheet are discussed. Accompanying the sheet 
development is the growth of a localised parallel electric field, one of the signatures of magnetic 
reconnection. Finally, the relevance of the results to a recent theory of turbulent reconnection is 
discussed. 
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I. INTRODUCTION 



The three-dimensional topological structure of many astrophysical plasmas, such as the 
solar corona, is known to be highly complex. In order to diagnose likely sites of energy re- 
lease and dynamic phenomena in such plasmas, where the magnetic Reynolds numbers are 
typically very large, it is crucial to understand at which locations strong current concentra- 
tions will form. These locations may be sites where singular currents are present under an 
ideal magnetohydrodynamic (MHD) evolution. In ideal MHD, the magnetic field is 'frozen 
into' the plasma, and plasma elements may move along field lines, but may not move across 
them. An equivalent statement of ideal evolution is that the magnetic flux through any 
material loop of plasma elements is conserved. 

Three-dimensional (3D) null points and separators (magnetic field lines which join two 
such nulls) might be sites of preferential current growth, in both the solar coronal and the 
Earth's magnetosphere, see Refs. 0,4] for reviews. 3D null points are thought to be present 
in abundance in the solar corona. A myriad of magnetic flux concentrations penetrate the 
solar surface, and it is predicted that for every 100 photospheric flux concentrations, there 
should be present between approximately 7 and 15 coronal null points^^. Furthermore, 
there is observational evidence that reconnection involving a 3D null point may be at work 
in some solar flares^ and solar eruptions*^. Closer to home, there has been a recent in situ 
observation^ by the Cluster spacecraft of a 3D magnetic null, which is proposed to play an 
important role in the observed signatures of reconnection within the Earth's magnetotail. 
In addition, current growth at 3D nulls has been observed in the laboratory^. 

While the relationship between reconnection at a separator (defined by a null-null line) 
and reconnection at a single 3D null point is not well understood, it is clear that the two 
should be linked in some way. Kinematic models anticipate that null points and separators 
are both locations where current singularities may form in ideal MHD^^. Moreover, it is 
also expected that 3D nulls^ 1 ^ and separators^ may each collapse to singularity in response 
to external motions. 

The linear field topology in the vicinity of a 3D magnetic null point may be examined by 
considering a Taylor expansion of B about the null; 

B = Mr, 

where the matrix Ai is the Jacobian of B evaluated at the null^ 1 ^. The eigenvalues of A4 



sum to zero since V ■ B = 0. The two eigenvectors corresponding to the eigenvalues with 
like-sign real parts define the fan (or S-) surface of the null, in which field lines approach 
(recede from) the null. The third eigenvector defines the orientation of the spine (or 7-) line, 
along which a single pair of field lines recede from (approach) the null [see Fig. [2(a)]. 

In Section [TT] we discuss the nature of the evolution of the magnetic field in the vicinity 
of a 3D null, and provide an example which demonstrates that certain evolutions of the 
null are prohibited in ideal MHD. In Sections II III and [TV] we present results of numerical 
simulations of a 3D null which is driven from the boundaries, and describe the qualitative 
and quantitative properties of the resulting current sheet. In Section [V] we observe that our 
models may point towards 3D nulls as a possible site where turbluent reconnection might 
take place. Finally in Section [VI] we give a summary. 

II. NON-IDEAL EVOLUTION AROUND 3D NULLS 
A. Ideal and non-ideal evolution 

The evolution of a magnetic field is said to be ideal if B can be viewed as being frozen 
into some ideal flow, i.e. there exists some w which satisfies 

E + wxB = -V$, (1) 

or equivalently 

<9B 

— = V x (w x B) (2) 

everywhere. In order for the evolution to be ideal, w should be continuous and smooth, such 
that it is equivalent to a real plasma flow. Examining the component of Eq. ([T]) parallel to B, 
it is clear that in configurations containing closed field lines, the constraint that § E ■ dl = 
must be satisfied in order to satisfy Eq. (pQ), so that $ is single- valued. However, even when 
no closed field lines are present, there are still configurations in which it may not be possible 
to find a smooth velocity w satisfying Eq. (TTJ). These configurations may contain isolated 
null points, or pairs of null points connected by separators^ 1 ^ 1 ^. We will concentrate in 
what follows on the case where isolated null points are present. 

Even when non-ideal terms are present in Ohm's law, it may still not be possible to find 
a smooth unique velocity w satisfying Eq. ([2]) if the non-ideal region is embedded in an 
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ideal region (i.e. is spatially localised in all three dimensions), as this imposes a 'boundary 
condition' on <3> within the (non-ideal) domain. If field lines link different points (say xi and 
x 2 ) on the boundary between the ideal and non-ideal region, then this imposes a constraint 
on allowable solutions for $ within the non-ideal region, since they must be consistent with 
$( Xl ) = $(x 2 ^. 

It has recently been claimed by Boozer— that the evolution of the magnetic field in 
the vicinity of an isolated 3D null point can always be viewed as ideal. However, we will 
demonstrate in the following section that in fact certain evolutions of a 3D null are prohibited 
under ideal MHD, and must therefore be facilitated by non-ideal processes. Furthermore, in 
subsequent sections we will show that such evolutions are a natural consequence of typical 
perturbations of a 3D null. 

The argument has been proposed^ that for generic 3D nulls 

(<9B \ 
—J (3) 

at the null point, where x is the position of the null, which can always be solved for w' 
so long as the matrix VB is invertible. For a 'generic' 3D null with det(VB) ^ this is 
always possible. However, Eq. (jHJ) is only valid at the null itself, with solution w' = dx^/dt. 
In general, V x (w' x B) ^ w' • VB in the vicinity of a given point. That is, it is not 
valid to discard the terms B(V • w') and (B ■ V)w' from the right-hand side of Eq. (j2J) 
(when expanded using the appropriate identity) when considering the behaviour of the flow 
around the null. In order to describe the evolution of the field around the null, we must 
consider the evolution in some finite (possibly infinitesimal) volume about the null. Under 
the assumption (J3j), the vicinity of the null moves like a 'solid body' — which seems to preclude 
any field evolution. Thus, the equation simply states that the field null remains, supposing 
no bifurcations occur — it shows that the null point cannot disappear, but does not describe 
the flux velocity in any finite volume around the null. This argument — that the velocity of 
the null and the value of the flux-transporting velocity at the null need not necessarily be 
the same — has been made regarding 2D X-points by Greene^. 

Furthermore, in order to show that the velocity w' is non-singular at the null, the as- 



sumption is made in Ref. 22j] that $ can be uniquely defined at the null, and that it may 
be approximated by a Taylor expansion about the null, which clearly presumes that <3> is 
a well-behaved function. As will be shown in what follows, in certain situations, it is not 



possible to find a function $ which is (or in particular whose gradient V$ is) well-behaved 
at the null point (smooth and continuous). This is demonstrated below by analysing the 
properties of the function $ with respect to any arbitrary boundary conditions away from 
the null point itself. By beginning at the null point and extrapolating $ outwards, it may 
indeed be possible to choose $ to be well-behaved at the null, but this $ must not be con- 
sistent with any physical (i.e. smooth) boundary conditions on physical quantities such as 
the plasma flow. 

In fact, it has been proven by Hornig and Schindler— that no smooth 'flux velocity' or 
velocity of 'ideal evolution' exists satisfying Eq. (j2J) if the ratios of the null point eigenvalues 
change in time. They have shown that under an ideal evolution, the eigenvectors (q a ) and 
eigenvalues (a) of the null evolve according to 

dq a da 

— = q a . Vw, - = Ca, (4) 

where C = — V • w| Xo is a constant. From the second equation it is straightforward to show 
that the ratio of any two of the eigenvalues must be constant in time. 

Note finally that the assumption that det(VB) ^ rules out the possibility of bifurcation 
of null points. While it is non-generic for det(VB) = to persist for a finite period of time, 
when multiple null points are created or annihilated in a bifurcation process, a higher order 
null point is present right at the point of bifurcation, and det(VB) passes through zero. 
Such bifurcation processes are naturally occurring, and clearly require a reconnection-like 
process to take place. 



B. Example 

We can gain significant insight by considering the kinematic problem. We consider an 
ideal situation, so that Ohm's law takes the form of Eq. ([1]). Uncurling Faraday's law gives 
E = — V$' — where B = V x A, and combining the two equations we have that 

~ dA 
V$ + wxB = — , (5) 
at 

where $ = $ — $'. We proceed as follows. We consider a given time-dependent magnetic 
field, and calculate the corresponding functions $ and Wj_. The component of w parallel to 
B is arbitrary with respect to the evolution of the magnetic flux. We choose the magnetic 
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field such that we have an analytical parameterisation of the field lines, obtained by solving 
= B(X(s)), where the parameter s runs along the field lines. Then taking the component 
of Eq. (jHJ) parallel to B we have 

V$"B = f-B 

$ = f%-Bds, 

where the spatial distance dl = |B|ds. Expressing B and ^ as functions of (Xo, s) allows us 
to perform the integration. This defines Vy$, while the constant of this integration, $ (X ), 
allows for different Vj_$. We now substitute in the inverse of the field line equations X (X) 
to find $(X). Finally, 

wi = * — 2 . (7) 

Consider now a null point, the ratios of whose eigenvalues change in time. Take the 
magnetic field to be 

B = (-2x-yf(t),y + xf(t),z). (8) 

A suitable choice for ^ is ^ = —f'/2 (0, 0, x 2 + y 2 ) (other appropriate choices lead to the 
same conclusions as below). The eigenvalues of the null point are 



A.= l , ^= - 1± f-^ . (9) 
with corresponding eigenvectors 

qi = (0,0,1) 

From the above it is clear that the spine (defined by (1, 0, 0) when / = 0) is given by 

z = , (l-A_)y + /x = 0. 

In addition, the equation of the fan plane (which lies at x = when / = 0) can be found by 
solving n ■ r = 0, n = qi x q + to give 
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Thus the spine and fan of the null close up on one another in time (being perpendicular 
when f{t) = 0), see Fig. QJb). Note that we require |/| < 3/2 to preserve the nature of the 
null point, otherwise is collapses to a non-generic null; we may take for example / = tanh(t). 
In order to simplify what follows, we define two new functions M and P which describe the 
locations of the spine and fan of the null: 

M(x, y,t) = {1 ')^ + fX , (10) 

P(x , yA = < i -y» + f x . (id 

1. Analytical field line equations 

Now, the first step in solving the equations is to find a representation of the field lines. 
Solving ^ = B(X(s)) gives 

x = (A + - l)C ie A+s + (A_ - l)C 2 e A - s . 

y = fC ie x + s + fC 2 e x ~ s (12) 
z = z e s 

where C\ and C 2 are constants. Clearly the equation for z is simple, but the other two are 
coupled in a more complicated way. In order to render the field line equations invertible, 
we choose to set s = on surfaces which move in time, tilting as the null point does. This 
is allowable since there is no linking between the integration of the field lines and the time 
derivative, that is, t (or /) is just a constant in the integration. We choose to set s = on 
surfaces defined by 

xif 

x = x 1 y , y = y - ^ - - (13) 

where X\ is a constant and y§ = yo(x,y,t) is the starting position of the field line footpoint 
on the tilting surface. Then 



We can see by comparison of equations ffTBj) and ffMj) that our 'boundaries' lie parallel to 
the fan plane, but with a shift of rkx\. 
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With C\ and C2 now denned by Eq. (fT4|) . Eqs. (fT2l) may be inverted to give 

z = zP- 1 ^-, (15) 

s = -^lnP. 

A_ 



2. Solving for <3? and physical quantities 

The first step in determining the solution is to solve for <3>, which is given by 

$ = f z (x 2 + y 2 ) ds (16) 



Jo 

where x, y, z are functions of (yo, zq,s). Upon substitution of (fl2]) . it is straightforward to 
carry out the integration in a symbolic computation package (such as Maple or Mathemat- 
ical since the integrand is simply a combination of exponentials in s. We then substitute 
( [151) into the result to obtain $(x, y, z). 

For any general choice of the function $0 (i-e. choice of initial conditions for the integra- 
tion, which may be viewed as playing the role of boundary conditions), we find that <3> is 
non-smooth at the fan plane, and therefore tends to infinity there. Thus the electric 
field E and flux velocity w will tend to infinity at the fan plane (see Fig. Ufa)). Examining 
the expression for $, it is apparent that the problematic terms are those in P~ l l x - (since 
— 1/A_ < —1/2 for — 1 < / < 1). In order to cancel out all of these terms, we must take (by 
inspection) 

s ° - -* ( a/- + a V + i_) »• < i7 » 

With this choice, $ is smooth and continuous everywhere. However, V$ is still non-smooth 
in the fan plane. Again examining $, we see that this results from terms in P^lnP (fii 
constant). These take the form (after simplification) 



$ = - ^fxtyozos 

Apx\ 



P- 1,X - MP~ X + /X - lnP" 1/A ". (18) 



A+ - A_ 

Now, due to the form of yo(x, y, z, t) and z Q (x,y, z,t) (Eq. fl5l) . it is impossible to remove 
the term in P m lnP through addition of any § (y , z ) without inserting some term in ei- 
ther z^lnz or M^lnM (/ij constant). This is simply equivalent to transferring the non- 
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FIG. 1: (a) Plots of w y against x for $0 = (black line) and for the case where <3?o is defined 
by Eq. (|17|) (grey line). The dashed line indicates the location of the fan plane, and we take 
/ = 0.2, x% = 1, y = 0.7, z = 0.5. (b) Orientation of the null spine and fan for t = (left) and t > 
(right). 

smoothness to the vicinity of the spine instead of the fan. Thus, there is no choice of 
boundary conditions which can possibly render the evolution 'ideal'. 

The example discussed above demonstrates that in an ideal plasma, it is not possible 
for a 3D null point to evolve in the way described, with the spine and fan opening/closing 
towards one another. Therefore, if such an evolution is to occur, non-ideal processes must 
become important. 



III. RESISTIVE MHD SIMULATIONS 



We now perform numerical simulations in the 3D resistive MHD model. The setup of 



yy Pontin and Galsgaard^. More details 



26IJ271 ]. We consider an isolated 3D null 



the simulations is very similar to that described 
of the numerical scheme may be found in Refs. 
point within our computational volume, which is driven from the boundary. We focus on 
the case where the null point is driven from its spine footpoints. The driving takes the form 
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FIG. 2: (a) (Colour online) Schematic of the 3D null point in the computational domain. The black 
arrows indicate the direction of the boundary driving, (b) Boundary driving flow at x = —Xl, the 
lower x-boundary, for Yi = Z\ = 3 and Ad = 80. 

of a shear. We begin initially with a potential null point; B = Bo(—2x,y, z), and with the 
density p — 1, and the internal energy e = 5/3/2 within the domain, so that we initially have 
an equilibrium. Here j3 is a constant which determines the plasma-/?, which is infinite at 
the null itself, but decreases away from it. We assume an ideal gas, with 7 = 5/3, and take 
B = 1 and (3 = 0.05 in each of the simulations. A stretched numerical grid is used to give 
higher resolution near the null (origin). Time units in the simulations are equivalent to the 
Alfven travel time across a unit length in a plasma of density p = 1 and uniform magnetic 
field of modulus 1. The resistivity is taken to be uniform, with its value being based upon 
the dimensions of the domain. Note that at t = 0, B is scale- free as it is linear, and thus, 
the actual value of r\ is somewhat arbitrary until we fix a physical length scale to associate 
with the size of our domain. 

At t = 0, the spine of the null point is coincident with the x-axis, and the fan plane 
with the x = plane [see Fig. [2(a)] . A driving velocity is then imposed on the (line-tied) 
x-boundaries, which advects the spine footpoints in opposite directions on the opposite 
boundaries (chosen to be in the ^/-direction without loss of generality). We choose an in- 
compressible (divergence-free) velocity field, defined by the stream function 




(19) 
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where y\ = yjY\ and z\ = z/Z\, the numerical domain has dimensions [±Xj, ±Y5, ±Z{\, and 
Ad describes the localisation of the driving patch. This drives the spine footpoints in the 
±y direction, but has return flows at larger radius (see Fig. G^b)). Note however that the 
field lines in the return flow regions never pass close to the fan (only field lines very close to 
the spine do). We have checked that there are no major differences for the evolution from 
the case of uni-directional driving (v = v y y). 

Two types of temporal variation for the driving are considered. In the first, the spine is 
driven until the resulting disturbance reaches the null, forming a current sheet (see below), 
then the driving switches off. It is smoothly ramped up and down to reduce sharp wavefront 
generation. In the other case, the driving is ramped up to a constant value and then held 
there (for as long as numerical artefacts allow). Specifically, we take either 

^) = t, »|(v) _1 ) ' °^^ r > ( 2 °) 

or 

K,(t) = t; tanh 2 (t/r), (21) 

where vq and r are constant. We begin by discussing the results of runs with the transient 
driving profile, as described by Eqs. f|T9|) and (1201) . 



A. Current evolution and plasma flow 

Unless otherwise stated, the following sections describe results of a simulation run with 
the transient driving time-dependence, and with parameters X\ = 0.5, Y\ = Zi = 3, Ad = 80, 
vq = 0.01, t = 1.8 and r] = 10 -4 . The resolution is 128 3 , and the grid spacing at the null is 
5x ~ 0.005 and 8y, Sz ~ 0.025 (the results have been checked at 256 3 resolution). 

As the boundary driving begins, current naturally develops near the driving regions. This 
disturbance propagates inwards towards the null at the local value of the Alfven speed. It 
concentrates at the null point itself, with the maximum of J = |J| growing sharply once 
the disturbance has reached the null, where the Alfven velocity vanishes. Figure [3] shows 
isosurfaces of J (at 50% of maximum). 

It is fruitful to examine the current evolution in a plane of constant z (plane of the 
shear). In such a plane, the current forms an 'S'-shape, and is suggestive of a collapse from 
an X-type configuration (spine and fan orthogonal) to a Y-type configuration in the shear 
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FIG. 3: Isosurfaces of J, at 50% of the maximum at that time, for times marked by diamonds in 
Fig. m (t = 1,2, 3, 5) 



plane (Fig. [5]). Note that once the driving switches off, the current begins to weaken again 
and spreads out in the fan plane, and the spine and fan relax back towards their initial 
perpendicular configuration (see Figs. [3]^d), E(d)). 

That we find a current concentration forming which is aligned not to the (global directions 
of the) fan or spine of the null (in contrast with simplified analytical self-similar— and 
incompressible^ solutions), but at some intermediate orientation, is entirely consistent with 
the laboratory observations of Bogdanov et al.—. In fact the angle that our sheet makes in the 
z = plane (with, say, the y-axis) is dependent on the strength of the driving (strength of the 
current), as well as the field structure of the 'background' null (as found in Ref. [jjj]) and the 
plasma parameters. We should point out that we consider the case where in the notation of 
Bogdanov et al. 7 < 1. They consider a magnetic field B = ((h+ h r )x, —{h — h r )y, —2h r z), 
and define 7 = h r /h. Thus, for 7 > 1 they have current parallel to the spine, which we 
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FIG. 4: Evolution of the maximum value of each component of J (J x dotted, J y dot-dashed, 
J z solid), as well as the integral of E» along the z-axis (dashed), and the time- variation of the 
boundary driving (triple-dot-dashed). The diamonds (A) mark the times at which the isosurfaces 
of J in Fig. [3] are plotted. 

expect to be resultant from rotational motions, and have very different current sheet and 
flow structures". 

Now examine the temporal evolution of the maximum value of each component of the 
current (see Fig. HI). It is clear that the current component that is enhanced significantly 
during the evolution is J z . This is the component which is parallel to the fan plane, and 
perpendicular to the plane of the shear (consistent with Ref. |25lj . and also with a collapse 
of the null's spine and fan towards each other (see, e.g. Ref. |19||)). 

Examining the plasma flow in the plane perpendicular to the shear, [z constant) we find 
that it develops a stagnation-point structure, which acts to close up the spine and fan (see 
Fig. [5]). It is worth noting that the Lorentz force acts to accelerate the plasma in this 
way, while the plasma pressure force acts against the acceleration of this flow (opposing 
the collapse). In the early stages after the current sheet has formed, the stagnation flow 
is still clearly in evidence, with plasma entering the current concentration along its 'long 
sides' and exiting at the short sides — looking very much like the standard 2D reconnection 
picture in this plane [Fig. El(b)]. As time goes on, the driving flow in the y-direction begins 
to dominate. 
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FIG. 5: (Colour online) Arrows show plasma flow, while shading shows |J| (scaled to the max- 
imum in each individual frame). Viewed in the z = plane, inner 1/4 of domain: [x, y] = 
[-0.12. .0.12,-0.7.-0.7] (y vertical). Images are at t = 1.6,2.4,3.0,5.0. 

B. Magnetic structure 

As the evolution proceeds and the current becomes strongly concentrated at the null, 
this is naturally where the magnetic field becomes stressed and distorted. Firstly, plotting 
representative field lines which thread the sheet and pass very close to the null, we see that 
in fact the topological nature of the null is preserved [Fig. [6](a)]. That is, underlying the 
Y-type structure of the current sheet there is still only a single null point present, with the 
angle between the spine and the fan drastically reduced (see below). 

It is interesting to consider the 3D structure of the current sheet, that is the nature of the 
quasi-discontinuity of the magnetic field at the current sheet. This is important since the 
common conception of a 'current sheet' is a 2D Green-Syrovatskii sheet with anti-parallel 



14 





-0.4 -0.2 0.0 ! 0.2 0.4 



0.6 



(c) 



-0.4 -0.2 



ill 


B 






>^v////' 


^^^^^ 



(d) 



FIG. 6: (Colour online) (a) Field lines traced from very close to the spine (black for negative x, 
grey (orange online) for positive), (b-d) The same, but for field lines traced from rings of larger 
radius around the spine and which graze the surface of the current sheet, viewed at different angles 
of rotation about the y-axis; (b) 7r/36, (c) 7t/6, (d) ir/2. This simulation run uses continuous 
driving, with parameters Xi = 0.5, Y\ = Z\ = 6, = 160, Vq = 0.04, r = 0.5 and r\ = 5 x 10~ 4 . 
Time of the images is t = 6. 



field on either side of a cut in the plane (see Refs. 3Q|]3l| ). Here, however, a current sheet 
localised in all three dimensions is present. 

In order to determine the structure of this 3D current sheet, we would like to know 
the nature of the jump (AB) in the magnetic field vector B across it, and how this varies 
throughout the sheet. Of course, because t] ^ 0, there is no true discontinuity in B, though 
the jump AB would occur over a shorter and shorter distance if r] were decreased (in an 
external domain of fixed size). This mismatching in B may be visualised by plotting field 
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lines which define approximately the boundaries of the current sheet, see Fig. [6](b-d). Along 
z = between the two ends of the current sheet, the black and grey (orange online) field 
lines (on opposite sides of the current concentration) are exactly anti-parallel [see Fig.E^d)]. 
Thus in this plane (only) we have something similar to the 2D picture (due to symmetry). 
However, for z ^ 0, the magnetic field vectors on either side of the current sheet are not 
exactly anti-parallel, and the black and grey (orange online) field lines cross at a finite angle 
(for small \y\). This angle decreases as \z\ increases (and as \y\ increases for z ^ 0), and 
thus the current modulus — proportional to AB across the sheet — falls off in y and z [see 
Fig. 0(d)], and is localised in all three dimensions [compare Fig. [3](c) and Fig. 0(c)]. 

C. Eigenvectors and eigenvalues 

Examining the evolution of the eigenvalues and eigenvectors of the null point in time 
provides an insight into the changing structure of the null. In order to simplify the discussion, 
we refer to the eigenvectors which lie along the x, y, and z axes at t = as the x-, y- and z- 
eigenvectors, respectively, and similarly for the eigenvalues. Consider first the eigenvectors. 
The orientation of the z-eigenvector is essentially unchanged, due to the symmetry of the 
driving. However, the x- and ^/-eigenvectors do change their orientations in time. In Fig. [71(a) 
the angle (6) between these two eigenvectors, or equivalently between the spine and fan, is 
plotted in time. In the case of continuous driving, the initial angle of 7r/2 quickly closes up 
to a basically constant value once the sheet forms. With the transient driving, the pattern 
is the same, except that the angle starts to grow again once the driving switches off (note 
that the minimum angle is smaller for the continually driven run plotted, since the driving 
Vq is 4 times larger). 

The change in time of the eigenvalues [see Fig. E(b)], and their ratios, is linked to that of 
the eigenvectors. The eigenvalues stay constant in time during the early evolution, but then 
begin to change, although they each follow the same pattern (i.e. they maintain their ratios, 
between approximately t = 1 and t = 2). By contrast, once the null point begins to collapse, 
there is a significant time dependence to the eigenvalues (t > 2), and evidently also to their 
ratios. This is clear evidence of non-ideal behaviour at the null, as described in Sec. HIl and 
demonstrates that the evolution modelled by the kinematic example given there is in fact 
a natural one. That is, the evolution prohibited under the restriction of ideal conditions is 
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FIG. 7: (a) Evolution of the angle (0) between the spine and fan (eigenvectors), for continual driving 
with parameters as in Fig. [6] (solid line) and the transient driving run (dashed), (b) Evolution of 
the null point eigenvalues (x-eigenvalue solid line, y dashed, z dotted), and the integrated parallel 
electric field along the z-axis (grey) for the transient driving run. The eigenvalues are normalised 
to their values at t = 0, for clarity. 

precisely that which occurs when the null point experiences a typical perturbation. 
D. Parallel electric field and reconnection 

Further indications that non-ideal processes are important at the null are given by the 
presence of a component of E parallel to B (E\\) which develops there. The presence of a 
current sheet at the null, together with a parallel electric field, is a strong indication that 
reconnection is taking place in the vicinity of the null (we actually expect that field lines 
change their connectivity everywhere within the volume defined by the current sheet (the 
'diffusion region'), not only at the null itself^i^, which is special only in that the footpoint 
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FIG. 8: Isosurface of En at t = 4.75 (time of its maximum), at 65% of maximum, transient driving 
run. 

mapping is discontinuous there). 



null. This quantity gives an exact measure of the rate of flux transfer across the fan surface 
when the non-ideal region is localised around the null point. The spatial distribution of E\\ 
within the domain when it reaches its temporal maximum is closely focused around the z- 
axis (see Fig. [8]). This is because J||B there by symmetry (note also that En is discontinuous 
at the plane z = by symmetry, since B z changes sign through this plane but J z (and thus 
E z ) is uni-directional — however, J E\\ ds is non-zero since ds also changes sign through 
z = 0). The field line which is coincident with the z-axis (by symmetry) thus provides the 
maximum value for J E\\ ds of any field line threading the current sheet — another reason to 
associate this quantity with the reconnection rate. 

The nature of the field line reconnection can be studied by integrating field lines threading 
'trace particles', which move in the ideal flow far from the current sheet, and which initially 
define coincident sets of field lines (see Fig. EJ). Once the current sheet forms, field lines 
traced from particles far out along the spine and far out along the fan are no longer coincident 
(i.e. they have been 'reconnected', near the null). Field lines traced from around the spine 
flip around the spine, while there is also clearly advection of magnetic flux across the fan 
surface. However, it seems that this occurs in the main part during the collapse of the null, 
while later in the simulations reconnection around/through the spine is dominant, since the 
boundary driving is across the spine. 




Endz provides a measure of the reconnection rate at the 
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FIG. 9: (Colour online) Field lines integrated from ideal trace particles, located initially far out 
along the spine (black) and the fan (grey, green and orange online). Images for (a) t = 0, (b) t = 3, 
(c) t = 6, and for continual driving with parameters as in Fig. [6J but vq = 0.03. 

IV. QUANTITATIVE CURRENT SHEET PROPERTIES 

In order to understand the nature of the current sheet that forms at the null point, its 
quantitative properties should be analysed. We focus here on two main aspects, firstly the 
scaling of the current sheet with the driving velocity, and secondly its behaviour at large 
times under continual driving. Crucial in this is whether the behaviour tends towards that 
of a Sweet-Parker- type current sheet, that is whether the sheet's length increases eventually 
to system-size, and whether the peak reconnection rate scales as a negative power of t], thus 
providing only slow reconnection for small t]. While the computational cost of a scaling 
study of Jmax with r\ for a fully 3D simulation such as ours is prohibitive, the investigations 
which we do perform provide an indication of the nature of the sheet. 

A. Scaling with vo 

Firstly we consider the scaling of the current sheet with the magnitude of the boundary 
driving velocity. In particular, we look at the maximum current density which is attained 
(which invariably occurs at the null), the maximum reconnection rate (calculated as the 
integrated parallel electric field along the fan field line coincident with the z-axis, as described 
previously), and also the current sheet dimensions L x , L y and L z (taken to be the full width 
at half maximum in each coordinate direction). This is done for the case of transient 
driving, with r fixed at a value of 1.8. One point which should be taken into account when 
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(a) (b) 

FIG. 10: (Colour online) Current modulus in the z = plane at time of current maximum, for (a) 
vq = 0.04, (b) vq = 0.001, scaled in each case to the individual maximum (20.4 in (a) and 0.44 in 
(b))- 

considering the measurements of the dimensions of the current sheet (particularly in the 
x- and y-directions) , is that the measurements do not necessarily mean exactly the same 
thing in the different simulation runs, in the sense that the current sheet morphology is 
not always qualitatively the same. Observe the difference between the current structures in 
Fig. [10], which shows cases with different driving strengths. When the driving in stronger, 
the current sheet is more strongly focused, and we measure a straight current sheet, which 
spans the spine and fan [Fig. flOlfa)]. However, when the driving is weaker, the region of 
high (above one half maximum) current density spreads along the fan plane, in an 'S' shape 
in and xy-plane [Fig. [3(c)]. Finally, when vq is decreased further still, we again have an 
approximately planar current sheet, but this time lying in the fan plane [Fig. fTDT b)]. Note 
that this changing morphology is also affected by the parameters i], (5 and 7, which we will 
discuss in a future paper, but which are held fixed here. 

We repeat the simulations with transient driving, varying vq, but fixing X\ = 0.5, Yi = 
Z\ = 3, Ad = 80, r = 1.8 and 77 = 5 x 10~ 4 . The results are shown in Fig. [TTJ First, we see 
that the peak current and peak reconnection rate both scale linearly with vq. The extent of 
the current sheet in x (L x ) increases linearly with v Q , which is a signature of the increased 
collapse of the null. By contrast, L y decreases with increased driving velocity. This is rather 
curious and seems counter-intuitive. It appears that the current sheet is more intense and 
more strongly focused at the null for stronger driving. In addition though, it is a result of the 
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peak reconnection rate (J* En), and the full width at half maximum of the current sheet in each 
coordinate direction (L x , L y , L z ). 



fact that an increasing amount of the current is able to be taken up in a straight spine-fan 
spanning current sheet, rather than spreading in the fan plane. Likewise the scaling of L z is 
also curious — the current sheet is more intense and strongly focused at the null for stronger 
driving. 

The above scaling analysis has also been performed for the case where the spine is dis- 
placed by the same amount each time, but at different rates. That is, as v is increased, r 
is decreased to compensate. In fact the scaling results are very similar to those above, but 
just with a slightly weaker dependence on vq. 
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B. Long-time growth under continual driving 

Now consider the case where, rather than imposing the boundary driving for only a 
limited time, the driving velocity is ramped up and then held constant. We take parameters 
as in Sec. IIII A[ but with Y\ — Z\ — 6 and run resolution 128 x 192 x 192, giving minimum 
8x ~ 0.005 and Sy, 5z ~ 0.025 again at the null. We focus on whether the current sheet 
continues to grow when it is continually driven, or whether it reaches a fixed length due 
to some self-limiting mechanism. As before, one of the major issues we run into between 
different simulation runs which use different parameters is in the geometry of the current 
sheet. For simplicity here we consider a case with sufficiently strong driving that the current 
sheet is approximately straight, spanning the spine and fan, as in Fig. [TOT a). 

There are many problems which make it hard to determine the time evolution of the 
current sheet length. During the initial stage of the sheet formation, it actually shrinks, as 
it intensifies, in the xy-plane (measuring the length as \JL\ + L 2 , see Fig. ITWa)). After this 
has occurred, the sheet then grows very slowly — differences over an Alfven time are on the 
order of the gridscale. There is no sign though of any evidence which points to a halting of 
this growth. Similarly, examining the evolution of L z , it seems to continue increasing for as 
long as we can run the simulations. Neither of the above is controlled by the dimensions of 
the numerical domain (see Fig. [TWb)). 

As to the value of the peak current, and the reconnection rate at the null, though the 
growth of these slows significantly as time progresses (see Fig. fl~2T c)). there is again no sign 
of them reaching a saturated value. The slowing of the current sheet growth (dimensions 
and modulus) is undoubtedly down to the diffusion and reconnection which occurs in the 
current sheet. 

In summary, it is difficult to be definitive that the current sheet will grow to system-size 
under continuous driving, due to the computational difficulties of running for a long time. 
However, there is no indication that the dimensions of the current sheet are limited by any 
self-regulating process. Rather, the dimensions are determined by the boundary conditions 
(i.e. the degree of shear imposed from the driving boundaries at a given time). Thus, if it 
were computationally possible to continue the shearing indefinitely, it appears that the sheet 
would continue to grow in size and intensity. If the system size is large and the resistivity 
is low, it is possible that the extended current sheet may break up into secondary islands, 
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(c) 

FIG. 12: Growth of the sheet in time, in (a) the xy-direction (y^L^ + Ly) and (b) the z-direction 
for three runs with different domain sizes (Yj = Zi = 6 solid line, Yj = Z\ = 3 dashed, Y\ = Z\ = 1.5 
dot-dashed), (c) Time evolution of the maximum value of each current component (J x solid line, 
J y dotted, J z dashed.) 

which lie beyond the scope of the present simulations. 

V. POSSIBILITY OF TURBULENT RECONNECTION 

A crucial aspect of any reconnection model which hopes to explain fast energy release 
is the scaling of the reconnection rate with the dissipation parameter. Mechanisms for 
turbulent reconnection have been put forward which predict reconnection rates which are 
completely independent of the resistivity^. Recent work by Eyink and Aluie^, who obtain 
conditions under which Alfven's "frozen flux" theorem may be violated in ideal plasmas, has 
placed rigourous constraints on models of turbulent reconnection. The breakdown of Alfven's 
theorem occurs over some length scale denned by the turbulence, which may be much larger 
than typical dissipative length scales. Eyink and Aluie demonstrate that in such a turbulent 
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FIG. 13: (Colour online) (a) Current and (b) vorticity profiles in the z = plane through the null, 
at the time of peak current, for parameters as in Sec. IIV Al with vq = 0.04. 

plasma, a necessary condition for such a breakdown to occur is that current and vortex 
sheets intersect one another—. This is a rather strong condition on the nonlinear dynamics 
underlying a turbulent MHD plasma. (Note that in order to correspond to any type of 
reconnection geometry, more than two current /vortex sheets should intersect, or equivalently 
they should each have multiple 'branches', as in Fig. [131 sa Y along the separatrices. This 
ensures a non-zero electric field at the intersection point, unlike in the simplified example of 
Eyink and Aluie.) 

One possible viewpoint of how such a situation might occur is that these current sheets 
and vortex sheets are generated by the turbulent mechanism itself. Alternatively, one might 
imagine another situation in which the result might be applicable is in a configuration where 
macro-scale current sheets and vortex sheets are present in the laminar solution, which may 
then be modulated by the presence of turbulence. The results discussed in the previous 
sections point towards 3D null points as sites where this might occur. 

Firstly, examining the vorticity (uj) profile in the simulations, we find that in fact a highly 
localised region of strong vorticity is indeed present. Furthermore, this region intersects with 
the region of high current density (see Fig. [T3l) and is focused at the null (and possibly spread 
along the fan surface as described above, depending on the choice of parameters). It also 
appears from the figure that in fact the vorticity profile forms in narrower layers than the 
current (in the xy-plane), one focused along the spine, with a much stronger layer in the 
fan (while the current sheet spans the spine and fan), so that rather than being completely 
coincident, the J and uj sheets really do 'intersect'. The rates at which J and uj fall off in z 
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away from the null are also very similar. 

Moreover, we have seen in Sec. Ill Bl that in order for a null point to evolve in certain 
ways, non-ideal processes are required. When the ideal system is supposed to evolve in such 
a way, a non-smooth velocity profile results, which typically shows up along either the spine 
or fan of the null, or both, depending on the boundary conditions. This non-smooth velocity 
corresponds to a singular vortex sheet. 

Finally, it is worth noting that in non-laminar magnetic fields, 3D nulls are expected to 
cluster together, creating 'bunches' of nulls^, thus providing the possibility that multiple 
coincident current and vortex sheets might be closely concentrated. 



VI. SUMMARY 



We have investigated the nature of the MHD evolution, and current sheet formation, at 
3D magnetic nulls. In complex 3D fields, isolated nulls are often considered less important 
sites of energetic phenomena than separator lines, due to the viewpoint that ideal MHD 
singularities in kinematic analyses are a result of the choice of boundary conditions. However, 
as demonstrated here (and proven in Ref. 2^])' certain evolutions of a 3D null are prohibited 
under ideal MHD, specifically those which correspond to a time dependence in the eigenvalue 
ratios of the null. We presented a particular example which demonstrates this, for the case 
where the angle between the spine and fan changes in time. The flow which advects the 
magnetic flux in the chosen example is shown to be non-smooth at the spine and/or fan for 
any choice of external boundary conditions. For typical boundary conditions (as required 
by say a line-tied boundary such as the solar photosphere, or at another null in the system) 
the flow will be singular. Thus non-ideal processes are always required to facilitate such an 
evolution. We presented the results of resistive MHD simulations which demonstrated that 
this evolution is a generic one. 

We went on to investigate the process of current sheet formation at such a 3D null in the 
simulations. The null was driven by shearing motions at the spine boundary, and a strong 
current concentration was found to result, focused at the null. The spine and fan of the 
null close up on one another, from their initial orthogonal configuration. Depending on the 
strength of the driving, the current sheet may be either spread along the fan of the null, or 
almost entirely contained in a spine-fan spanning sheet (for stronger driving). The structure 
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of the sheet is of exactly anti-parallel field lines in the shear plane at the null, with the 
intensity of J falling off in the perpendicular direction due to the linearly increasing field 
component in that direction, which is continuous across the current sheet. Repeating our 
simulations but driving at the fan footpoints instead of the spine, a very similar evolution 
is observed. The current is still inclined to spread along the fan rather than the spine for 
weaker driving, consistent with relaxation simulation results^. This is natural due to the 
shear driving, and the disparity in the structures of the spine and fan. The spine is a single 
line to which field lines converge, and the natural way to form a current sheet there is to 
have those field lines spiral tightly around the spine line. This field structure however is 
synonymous with a current directed parallel to the spine, and is known to be induced by 
rotational motions, rather than shearing ones^. 

In addition to the current sheet at the null, indications that non-ideal processes and 
reconnection should take place there are given by the presence of a localised parallel electric 
field. The maxima of this parallel electric field occur along the axis perpendicular to the 
shear (z). The integral of E\\ along the field line which is coincident with this axis (by 
symmetry) gives the reconnection rate^, which grows to a peak value in time before falling 
off after the driving is switched off. Field lines are reconnected both across the fan and 
around the spine as the null collapses. 

As to the quantitative properties of the current sheet, its intensity and dimensions are 
linearly dependent on the modulus of the boundary driving. In addition, under continual 
driving, we find no indications which suggest that the sheet does not continue to grow (in 
intensity and length) in time. That is, its length does not appear to be controlled by any 
self-regulating mechanism. Rather, its dimensions at a given time are dependent on the 
boundary conditions (the degree of shearing). 

Finally, with respect to a recent theorem of Eyink and Aluie^, our results suggest 3D 
nulls as a possible site of turbulent reconnection. Their theorem states that turbulent 
reconnection may occur at a rate independent of the resistivity only where current sheets 
and vortex sheets intersect one another. We find that strong vorticity concentrations are in 
fact present in our simulations, and are localised at the null and clearly intersect the current 
sheet. Moreover, in the ideal limit, our kinematic model demonstrates that as the spine 
and fan collapse towards one another, the vorticity at the fan (or spine, or both) must be 
singular (due to the non-smooth transport velocity). 
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